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Abstract. When two Markov operators commute, it suggests that we can 
couple two copies of one of the corresponding processes. We explicitly con- 
struct a number of couplings of this type for a commuting family of Markov 
processes on the set of conjugacy classes of the unitary group, using a dynam- 
ical rule inspired by the RSK algorithm. Our motivation for doing this is to 
develop a parallel programme, on the circle, to some recently discovered con- 
nections in random matrix theory between reflected and conditioned systems 
of particles on the line. One of the Markov chains we consider gives rise to 
a family of Gibbs measures on 'bead configurations' on the infinite cylinder. 
We show that these measures have determinantal structure and compute the 
corresponding space-time correlation kernel. 



1. Introduction 

When two Markov operators commute, it suggests that we can couple two copies 
of one of the corresponding processes. Such couplings have been described in |10j 
in a general context. In this paper, we explicitly construct a number of couplings of 
this type for a commuting family of Markov processes on the set of conjugacy classes 
of the unitary group, using a dynamical rule inspired by the Robinson-Schensted- 
Knuth (RSK) algorithm. Our motivation for doing this is to develop a parallel 
programme, on the circle, to some recently discovered connections in random matrix 
theory between reflected and conditioned systems of particles on the line (see, for 
example, [2j [4], I19L I32L I40j ). The RSK algorithm is a combinatorial device which 
plays an important role in the representation theory of GL(n,C) and lies at the 
heart of these developments. We refer the reader to |16] for more background on 
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the combinatorics. One of the Markov chains we consider gives rise to a family of 
Gibbs measures on 'bead configurations' on the infinite cylinder. This is related 
to recent work [51 [6j 124] on planar and toroidal models. We will show that these 
measures have determinantal structure and compute the corresponding space-time 
correlation kernel. 

We start with some motivation, and a flavour of some of the main results in 
this paper. The following construction is closely related to the RSK algorithm [321 
EMS M\- Let Si = Z and, for n > 2, 

E n = {x e Z" : x 1 < ■■■ < x„}. 

The reader should think of an element x € E n as a configuration of n particles 
on the integers located at positions x\ < • • • < x n . We say that a pair (x,y) € 
E m x E m+ i are interlaced, and write x ^ y, if yj < Xj < yj+i for all j < m. A 
(discrete) Gelfand-Tsetlin pattern of depth n is a collection (x 1 ^ 2 , ...,x n ) such 
that x m £ E m for m < n and x m < x m+1 for 1 < m < n. 

Fix n > 2 and denote by GT n the set of Gelfand-Tsetlin patterns of depth n. Let 
ioi, W2, ■ ■ ■ be a sequence of independent random variables, each chosen according to 
the uniform distribution on {1, 2, . . . , n}. Using these, we will construct a Markov 
chain (X(k), k > 0) with state space GT„, which evolves according to 

X(k + l) = g(X(k),w k+1 ), k>0 

where g : GT n x {1, . . . ,n} — > GT„ is defined recursively as follows. Fix 77i < n 
and let [x, y) £ £" m x E m+ \ such that x ^ y. Let cc' S £" m such that, for some 
j < m, x[ — Xi + Sij. The reader should have in mind m particles located at 
positions x± < ■ ■ ■ < x m , interlaced with another set of m + 1 particles located 
at positions y\ < ■■■ < y m +i- The first configuration x is updated by moving 
the particle at position Xj one step to the right, that is, to position Xj + 1, giving 
a new configuration x'. This can be used to obtain an update y' to the second 
configuration y, obtained by moving the first available particle, strictly to the right 
of position yj, which can be moved without breaking the interlacing constraint, 
so that x' ^ y' . Such a particle is guaranteed to exist because the interlacing 
constraint cannot be broken by the rightmost particle. In other words, the updated 
configuration is given by y[ = y-i + 5ik where k = inf{Z > j : yi + 1 < x{\. Let 
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us write y' = <j)(x,y,x') where <f> is defined on an appropriate domain. Now, given 
rn < n and a pattern (x 1 , . . . , x n ) G GT„ we define a new pattern 

(y 1 ,...,y n ) = g((x 1 ,...,x n ),m) 

as follows. First we set y J = x J for j < m. Then we obtain y m from the configura- 
tion x m by moving the first available particle, starting from the particle at position 
x™, which can be moved one step to the right without breaking the interlacing 
constraint. Finally, we define, for m < I < n, y l+1 — (j)(x , x l+1 , y). 

The Markov chain X has remarkable properties. For example, if we start from 
the initial pattern 

(1) X(0) - ((0), (-1, 0), (-2, -1, 0), . . . , (-n + 1, . . . , 0)), 

then (X n (k), k > 0) is a Markov chain (with respect to its own filtration) with state 
space E n and transition probabilities given by 



(2) P n (x,y) 



1HV1 x / y 

n h(x) ' " 

otherwise 



where h is the Vandermonde function h(x) — Y[i<j< n ( x 3 ~ Xi ) an< ^ ^ ne notation 
x y y means that the configuration y can be obtained from x by moving one particle 
one step to the right. The fact that P n is a Markov transition kernel follows from 
the fact (see, for example, |25| ) that h is positive on E n and satisfies 

±-"£h(y) = h(x). 

More generally (see, for example, [301 131p : 

Proposition 1.1. Given x £ E n> if X(0) is chosen uniformly at random from 
the set of patterns (x , . . . , x n ) with x n — x, then (X n (fe), k > 0) is a Markov chain 
started from x with transition probabilities given by P n . Moreover, for each T > 0. 
the conditional law of X(T), given (X n (k),T > k > 0), is uniformly distributed on 
the set of patterns (x 1 , . . . , x n ) with x n = X n (T). 

The connection with the RSK algorithm is the following: if we start from the 
initial pattern |T]), then the integer partition 

(X^klX^k) + l,x:_ 2 (k) + 2,..., X?{k) +n - 1) 
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is precisely the shape of the tableau obtained when one applies the RSK algorithm, 
with column insertion, to the word W1W2 ■ ■ ■ Wk- We refer to [30] for more details. 
This construction clearly has a nested structure. If we consider the evolution of 
the last two rows X™ -1 and X n we see that we can construct a Markov chain with 
transition probabilities P n from a Markov chain with transition probabilities P n -i 
plus a 'little bit of extra randomness'. 

In the above construction we are thinking in terms of particles moving on a 
line. In random matrix theory, there are often strong parallels between natural 
measures on configurations of particles (or 'eigenvalues') on the line, and configura- 
tions of particles on the circle. It is therefore natural to ask if there is an analogue 
of the above construction for configurations of particles on the circle. The notion 
of interlacing carries over in the obvious way. However, in this setting, interlaced 
configurations should have the same number of particles, and the analogue of a 
Gelfand-Tsetlin pattern could potentially be an infinite object. Despite this fun- 
damental difference between the two settings, there is indeed a natural analogue 
of the above 'RSK dynamics' and a natural analogue of Proposition ll.il Consider 
the discrete circle with N positions which we label {0, 1, . . . , N — 1} in the anti- 
clockwise direction. The analogue of the Markov chain with transition matrix P n 
is a Markov chain on the set of configurations of n indistinguishable particles 
on the discrete circle with transition probabilities given by 



where, similarly as before, the notation x /* y means that the configuration y can 
be obtained from x by moving one particle one step anti-clockwise, and the function 
A is again a Vandermonde function defined, for a configuration x which consists of 
a particles located at positions ki , . . . , k n , by 



(3) 




otherwise 



(4) 



A(x)= J] \e 




i<j<n 



The constant c is chosen so that Q(x, •) is a probability distribution; the fact that 
c can be chosen independently of x follows from the fact (see, for example, [25 ) 
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that A is a positive eigenvector of the matrix l x yyi that is, 

£ A(y) = XA(x), 

for some A > 0. A formula for the eigenvalue A can be found in [25) . 

In this setting we will say that a pair of configurations (x, y) are interlaced, 
and write x -< y as before, if there is a labelling of the particles such that x consists 
of a particles located at positions k\, . . . ,k n , y consists of a particles located at 
positions l\, . . . ,l n , kj < lj < fey+i for j < n and either k n < l n < k\ + N or 
k n < l n + N < ki + N. 

Let (X(k),k > 0) be a Markov chain with state space and transition 
matrix Q. On the same probability space, without using any extra randomness, 
we can construct a second process (Y(k),k > 0), also taking values in , such 
that X(k) ^ Y(k) for all k. This is given as follows. For each k > 0, given X(k), 
Y(k) and X(k + 1) we obtain the configuration Y(k + 1) from Y{k) by moving a 
particle one step anti-clockwise; this particle is chosen as follows. The transition 
from X{k) to X(k + 1) involves one particle moving one step anti-clockwise; starting 
at the position of this particle, choose the first particle in the configuration Y(k) 
that we come to, in an anti-clockwise direction, which can be moved by one step 
anti-clockwise without breaking the interlacing constraint. 

The function A defined by (J4J) is also a positive (left and right) eigenvector of 
the matrix l x -<y This follows, for example, from the discussion given at the end of 
Section 3. In particular, 

^A(x) = 7 A(y) 

for some 7 > and we can define a Markov kernel on by 

1 AO) 

(5) M(y,x) = -^l x ^ y . 

The analogue of Proposition 11.11 in this setting is the following: 

Proposition 1.2. IfY(0) = y and X(0) is chosen at random according to 
the distribution M(y, •), then (Y(k),k > 0) is a Markov chain started at y with 
transition matrix Q. Moreover, for each T > 0, the conditional law of X(T), given 
(Y(k), T > k > 0), is given by M(Y(T), ■). 
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In the sequel we will present a number of variations of this result, firstly involv- 
ing random walks with jumps in a continuous state space and secondly involving 
Brownian motion. Proposition 11.21 follows from results presented in Section 3 (see 
discussion towards the end of that section) . 

We will also study continuous analogues of the Markov chain with transition 
matrix M. These Markov chains also arise naturally in the context of a certain 
random walk on the unitary group which is obtained by taking products of certain 
random (complex) reflections, as studied for example in [15|, 135] . This is described 
in Section 2 and taken as a starting point for the exposition that follows. The main 
point is that these Markov chains commute with each other and with the Dirichlct 
Laplacian on the set of conjugacy classes of the unitary group. In Sections 3 and 4, 
we present couplings which realise these commutation relations, first between inter- 
laced random walks on the circle and later between interlaced Brownian motions. 
These couplings are precisely the variations of Proposition 11.21 mentioned above. 
Actually there are two natural couplings for the random walks and these corre- 
spond to dynamical rules inspired by the RSK algorithm with row, and column 
insertion, respectively. The couplings between interlaced Brownian motions can be 
thought of as a limiting case where the two types of coupling become equivalent. 
In Section 5, we consider a family of Markov chains which can be thought of as 
perurbations of a continuous analogue of the Markov chain with transition matrix 
given by ([5]). These give rise to a natural family of Gibbs measures on 'bead config- 
urations' on the infinite cylinder. This is a cylindrical analogue of the planar bead 
model studied in [5] and, in one special case (the 'unperturbed' case), can also be 
regarded as a cylindrical analogue of some natural measures on Gelfand-Tsetlin 
patterns related to 'GUE minors' [2], 1231 1331 114] (see also [7] for extensions to the 
other classical complex Lie algebras). We show that these measures have determi- 
nantal structure by first writing the restrictions of these measures to cylinder sets 
as products of determinants and then following the methodology of Johansson (see, 
for example, [22 ) to compute the space-time correlation functions. 

2. Markov processes on the conjugacy classes of the unitary group 

Consider the group U(n) of n x n unitary matrices, and denote by C n the 
set of conjugacy classes in U(n). Each element of C n can be identified with an 
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unlabclled configuration of n points (eigenvalues) on the unit circle, which in turn 
can be identified with the Euclidean set D n = {0 S W 1 : < 6\ < ■ ■ ■ < 6 n < 2tt}. 
Denote by dx the image of Lebesgue measure under the latter identification, and 
by ji the probability measure on C„ induced from Haar measure on U(n). Then 
li(dx) = (2ir)~ n A(x) 2 dx, where A(x) is defined, for x = {e ldl , . . . , e l9 ™}, by 



Kl<m<n 



(6) A(x) = n 

l<Km< 

The irreducible characters x\ of U(n) are indexed by the set 

n n = {\EZ n : Ai > A 2 > • • • > A„}. 

We assume that these are normalised, that is, J |xa| 2 ^ = 1 for all A £ il n . 

Let x S C n and consider the random walk on U (n) which is constructed by 
multiplying together independent, randomly chosen elements from the conjugacy 
class x. (By 'a randomly chosen element from the conjugacy class a;' we mean a 
random element of the conjugacy class x which admits a representation of the form 
MD X M*, where D x is an arbitrary element of the conjugacy class x and M is a 
Haar-distributed, randomly chosen element of U(n).) The corresponding Markov 
kernel p x (y,dz) can be interpreted as the law of D x MD y M* , where D x and D y 
are arbitrary elements of the conjugacy classes x and y, respectively, and M is a 
Haar-distributed randomly chosen element of U(n). Write p x f = J p{-, dz)f(z) for 
/ € L 2 (C„,^). Then, for each A e Q n , 

/-n, Xx(x) 
(<) PxX\ = —, — XX, 

d\ 

where d\ is the dimension of the representation corresponding to A (see, for exam- 
ple, |13[ Proposition 6.5.2]). Moreover, fj, is an invariant measure for the Markov 
kernel p x , that is, 

(8) J n(dy)p x (y, ■) = fi. 

The Markov kernels {p x , x € C n } are the extreme points in the convex set M. 
of all Markov kernels on L2(C„, fi) with the irreducible characters as eigenfunctions. 
(See, for example, PQ).) All of the Markov kernels in Ai have /i as an invariant mea- 
sure and, as operators on Li(C n ,fi)^ they commute. They correspond to random 
walks on U (n) such that the law of the increments is invariant under conjugation. 
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Such random walks on U(n), and other classical compact groups, have been studied 
extensively in the literature (see, for example, [3 [H \38\ 1341 135j), 

The case of interest in this paper is the random walk obtained by taking prod- 
ucts of certain random (complex) reflections in U(n) (see, for example, I15|, I35j ). 
More precisely, we take x = {e lr , 1, 1, . . . , 1} in the above kernel, where r e (0, 2w). 
Let us write 

(9) p r := p x 

for this case. A concrete description of this kernel can be given as follows (see |15j 
for details). For a, b 6 D n , write a < b if 

a\ < b\ < a% < ■ ■ ■ < a n < b n . 

For y — {e lai , . . . , e M "} and z — {e lbl , . . . , e lbn }, where a, b <E D n , write y < r z if 
either a <b and ^2j(bj — a,j) — r, or b < a and ^2j(bj — aj) + 2tt = r. The measure 
p r (y, ■) is supported on the set 

F r {y) = {z e C„ : y < r z}. 

This set can be identified with the disjoint union of a pair of (n — l)-dimensional 
Euclidean sets 

{b e D n : a<b, - aj) = r} U {6 e D n : a^b, - = r - 27r}, 

J 3 

each of which can be endowed with (n — l)-dimensional Lebesgue measure giving a 
natural measure on their union. The measure obtained on F r (y) via this identifi- 
cation can be extended to a measure v r {y, •) on C n by setting v r (y, C n \F r (y)) = 0. 
The following identity can be deduced from [151 Lemma 2]. 

Proposition 2.1. 

(10) Pr(y,dz) = — ^-i/ r (y,dz), 

lr A(y) 

where j r = |1- e lr \( n ~^ / (n - 

In Section 5, we will present a determinantal formula for q r v r (y,dz)dr, 
where q > is a parameter, and use this to give an alternative proof of Proposi- 
tion o 
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Another operator which will play a role in this paper is the Dirichlct Laplacian 
on the (closed) alcove 

(11) A n = {9 e M" : 0i < 2 < • • • < On < 0i + 2tt}. 

Let (Qt) denote the transition semigroup of a standard Brownian motion condi- 
tioned never to exit A n . This is a Doob transform of the Brownian motion which 
is killed when it exits A n , via the positive eigenfunction 

(12) h(9)= J] |e i9i -e i( H 

l</<m<n 

We can identify C n with exp(z^4 n )/C„, where £ n denotes the group of cyclic per- 
mutations which acts on A n by permuting coordinates. It is known [3] that the 
eigenfunctions of the induced semigroup (Qt) on Li(C n ,[i) are given by the irre- 
ducible characters {x\, A S £l n }, which implies that Q t G M, for each t > 0. Note 
that the corresponding process on C n can be thought of as n standard Brownian 
motions on the circle conditioned never to collide. 

We will also consider discrete analogues of the above processes. Set 

A% = (NA n /2ir)nZ n , = eM^iA n (N)/N)/£ n , fi^{Ae Q n : X 1 < N-l}. 

Much of the above discussion can be replicated in this setting, but for our purposes 
it suffices to make the following remark. Think of as the set of configurations 
of n particles at distinct locations on the discrete circle with N positions. Consider 
the random walk in where at each step a particle is chosen at random and 
moved one position anti-clockwise if that position is vacant; if it is not vacant the 
process is killed. It is known that the restriction of the function A to is the 
Perron- Frobenius eigenfunction for this sub-Markov chain. In fact, a complete set 
of eigenfunctions (with respect to the measure A(x) 2 ) is given by the restrictions 
of the characters {x\, A € f2^} to C% ■ (See, for example, |25j.) As we shall see 
later, the discrete analogues of the v r (thought of as operators) commute with the 
transition kernel of this killed random walk and therefore share these eigenfunctions. 

3. Couplings of interlaced random walks 

The Markov chain on C n with transition probabilities p r , defined by Q, can be 
lifted to a Markov chain on A n , which is better suited to the constructions of this 
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section. To make this precise let us say x and x' belonging to A n are r-interlaced, 
for some r € (0, 2n), if 

n 

(13) x- e [xi, Xi+i] for i = 1,2, . . . ,n and /J(x^ — £i) = r, 

t=i 

when we adopt the convention that Xn+i = x\ + 27r. In this case we will write 
x < r x' . Define 7r : A n — » C„ by 7r(x) = {e lXl , . . . , e lXn }. Denote by l r (x, dx') the 
(n — l)-dimensional Lebesgue measure on the set 

G r (x) = {x 1 € A„ : x ;< r x'}. 

Clearly the restriction of 7r to G r (x) is injective, with 7r(G r (x)) = F r (7r(x)). More- 
over, for measurable B C A n , l r (x, B) = v r (ir(x), Tr(BnG r (x))). Thus, if we define, 
for measurable B C A„ , 

(14) q r (x, B) = Pr(Tr(x), ir(B fl G r (x))) 
then, by Proposition ^. II 

q r {x,B) = / p r (-n(x),dz) 

Jir(BnG r (x)) 

Jir{BnG T {x)) lr A(tt(x)) 

- J Blr A(n(x)) lr[X ' aXh 

and hence, 

(15) K^-i^M**-). 

7r A(7T(x)) 

We will refer to a Markov chain with values in A„ and transition probabilities 
q r as an r-interlacing random walk. As far as we are aware, such processes have 
not previously appeared in the literature. Note that since p r p s = PsPr we have 
<7r<7;s — 9s9r or i equivalently, l r l s — l s l r . 

The goal of this section is to construct, for given r, s S (0, 2tt), two different 
Markovian couplings (X(k),Y(k);k > 0) of a pair of r-interlacing random walks 
on A n , having the property that X(k) and Y{k) are s- interlaced for each n and 
moreover, for each / > 0, the trajectory (Y(k);0 < k < I) will be a deterministic 
function of Y(0) together with the trajectory (X(fc); < k < I). 
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The existence of such couplings is suggested by the commutation relation q r q s = 
q s q r , equivalently l r l s = l s l r . For any u,v € A n consider the two sets t u ,v = {x S 
A n : u < s x ^ r v} and t' uv = {y £ A n : u ^ r y ^ s v}. If either is non-empty, 
then they both are, and in this case they are (n— l)-dimensional polygons, and the 
relation l r l s = l s l r can be interpreted as saying these two polygons have the same 
(n — l)-dimensional volume. In fact the two polygons are congruent. Define 

(16) y = 4>u,v{x) 

via 

(17) yi = VDm{u i+ i,Vi) + max(uj, Vi-i) - x t . 

It is easy to see that <j> UiV is an isometry from r UtV to t' u v using the facts that 
yi e [max(tif,i;t_i),min(tti+i,i;t)] if and only if n € [max(tti,tii_i),min(ii<+i,Ui)], 
and 

n n n 

^2vi = y^(min(Mj + i,^) + max(uj,Uj_x) - Xi) = J^(w 4 + «i - Xi). 

i—l i—1 i—1 

Proposition 3.1. Let (X(k); k > 0) be an r -interlacing random walk, starting 
from X(0) having the distribution q s (y,dx) for some given y € A n , where q s is 
defined by Jj^[ ). Let the process (Y(k); k > 0) be given by Y(0) = y and 

Y(k + 1) - 4> Y (k),x(k + i){X{k)), for k > 0, 

where <f> is defined by il6}) . Then (Y{k);k > 0) is distributed as an r -interlacing 
random walk starting from y. 

Proof. We prove by induction on m that the law of (Y(l), . . . Y(m), X(m)) is 
given by q r (y, dy(l)) . . . q r (y(m— 1), dy(m))q s (y(m), dx(m)). Suppose this holds for 
some m. Then, since Y(l) . . . Y(m) are measurable with respect to ^(0), X(l), . . . X{m) 
the joint law of (Y(l), . . . Y(m), X(m), X(m + 1)) is given by 

Qr(y,dy(l)) ■ . . q r (y(m - l),dy(m))q s (y(m),dx(m))q r (x(m),dx(m + 1)). 

Equivalently we may say that the law of (Y(l), . . . Y(m),X(m + 1)) is 



Qr{y,dy(l)) . . . q r (y(m - l),dy{m))(q s q r )(y(m) 7 dx(m + 1)) 
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and that the conditional law of X(m) given the same variables is uniform on 
T Y(m)x(m+i)- From the measure preserving properties of the maps (j> u ,v, it fol- 
lows that, conditionally on (Y(l), . . . Y(m), X(m + 1)), Y(m + 1) is distributed 
uniformly on Tw TO \ x(m+i)> an< ^ ^he inductive hypothesis for m + 1 follows from 
this and the commutation relation q s q r = q r q s . □ 



The dynamics of the coupled processes (X(k), Y(k); k > 0) are illustrated in the 
following two diagrams, in which interlacing configurations x = X{k) and y = Y(k) 
are shown together with updated configurations x' — X(k + 1) and y' = Y(k + 1). 
It is natural to think of these as particle positions on ( a portion of ) the circle. For 
simplicity we consider an example where x' and x differ only in the i-th co-ordinate. 
Vi-i Xi-i yi Xi y l+ i x i+ i 

f t \\ f — r~ 

-0; 9; + & #^ 

Vi-l x i-l Vi x i Vi+l x i+l 

The configuration y' is of course determined by y, x and x' together. The 
simplest possiblity is shown above, in this case ith y particle advances by the same 
amount as the zth x-particle. However should the zth x-particle advance beyond 
yi+i then it it pushes the (i + l)th y particle along, whilst the increment in the 
position of the ith y particle is limited to y%-\-\ — X{. 



Vi-x yi Xi y i+ i x l+1 




The proof of Propostion l3.1l made use of only the measure preserving properties 
of the family of maps <p uv , and consequently we can replace it by another family 
of maps having the same property and obtain a different coupling of the same 
processes. The pushing interaction in the coupling constructed above is also seen 
in the dynamics induced on Gelfand-Tsetlin patterns by the RSK correpondence. 
There exists a variant of the RSK algorithm (with column insertion replacing the 
more common row insertion), in which pushing is replaced by blocking. We will 
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next describe a family ip uv of measure-preserving maps that lead to a coupling with 
such a blocking interaction. 

We recall first a version of the standard Skorohod lemma for periodic sequences. 

Lemma 3.2 (Skorohod). Suppose that (zi\i G Z) is n-periodic and satisfies 
Y^i=i z i < 0- Then there exists a unique pair of n-periodic sequences (r,;i G Z) 
and (k; G Z) such that 

r i+ i = ri + Zi + k+i for all i G Z, 

mi/i i/ie additional properties ri >0,h> 0, and k > => = for all i G Z. 

A configuration (x\^X2 ...x n ) of n points on the circle will be implicitly ex- 
tended to a sequence (xi\i G Z) satisfying = Xi + 2n. 

PROPOSITION 3.3. Suppose that u, v, and x are three configurations on the 
circle with u < s x -< r v, where s > r. Define an n-periodic sequence via 

z l = v l - x l - x i+1 + u l+1 for i G Z, 

and let (r,l) be the associated solution to the Skorohod problem. Set yi — Xi — l i} 
then (j/i, . . . , y n ) is configuration on the circle such that u < r y < s v. 

PROOF. Notice that Y2i=i Zi = r — s < so the appeal to the Skorohod con- 
struction is legitimate. 

Summing the Skorohod equation gives Y2i=i Zi = ~ 127=1^- Consequently, 
using the definition of y, 

n n n 

y^(yi - «i) = y^( x i - h - u i) = s - y^ji = r - 

i—l i—l i—l 

Similarly ELi(^ _ Vi) = s - 

To conclude we will verify the inequalities 

max(u J ,i; J _i) <y i <x l < min(ti i+ i, Vi). 

It is easily checked that z~_ 1 < min(a; i — Ui, Xi — Vi-i) . The solution of the Skorohod 
problem satisfies < h < z~_ 1 . Together with the definition of j/j, this gives the 
desired inequalities. □ 
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By virtue of the preceeding result we may define ip u ,v ■ i~ uv — ► t' uv via ip U v(x) = 

y- 

Proposition 3.4. ip uv is a measure-preserving bisection. 

Proof. For x e A n , extended as before to a sequence (xf,i € Z, let x^ be 
denned by x\ = —x-i for igZ. Observe that if x -< r y, then yt p or Uj Vj ^ 

and y as in the previous proposition we will show that the application 

(u,x,v) i ^ (u t ,y t ,u t ) 

is an involution, and this implies in particular that ip u , v is invertible with art = 
"4> v t u i(y^)- To this end first note that ~< s ~< r u\ and hence it is meaningful to 
apply VVut to y t . We must proof that the resulting conguration x say, is equal to 
art. 

We have ar, = y\ — U, where (f, Z) solves the Skorohod problem with data 

Zi=u\-y\-y\ +l +v\ +1 . 

Since ar, — y\ — h = —X-i + l-i — h verifying that x = art boils down to checking 
that U = l-i. For this it is suffices, by the uniqueness property of solutions to 
the Skorohod problem, to confirm that r[ — r_j and l\ = l-i solve the Skorohod 
problem with data z. Now r' and /' arc non-negative and satisfy l[ > ==>■ r[ — 
because r and I have these properties. We just have to compute 

r'i + Zi + = r- t + u\ - y\ - y\ +1 + v\ +1 + 

= r-i - u-i + y-i + y-( i+ i) - t>-(»+i) + 

= r-i - u-i + x-i - l- t + ar_ (i+1) - l-a+i) - v-( i+ i) + 

= r-i - z-(i+i) - l-i 

= r -(i+i) = r'i+i- 

This proves the involution property. 

It remains to verify the measure-preserving property. The construction of tp uv 
is such that it is evident that it is a piecewise linear mapping, and that its Jacobian 
is almost everywhere integer valued. Since the same applies to the inverse map 
constructed from V„t u t we conclude the Jacobian is almost everywhere ±1 valued. 

□ 
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Finally the following propostion follows by exactly the same argument as Propo- 
sition 13.11 



Proposition 3.5. Suppose that s > r. Let (X(k);k > 0) be an r --interlacing 
random walk, starting from X(0) having the distribution q s (y,dx) for some given 
y G A n . Let the process (Y(k); k > 0) be given by Y(0) — y and 

Y(k + 1) = i> Y (k),x(k+i) {X(k)). for k>0. 

Then (Y(k)]k > 0) is distributed as an r -interlacing random walk starting from y. 

Let us illustrate this coupling in a similar fashion to before. 

yi x l y i+ i x i+ i y i+2 x l+2 



-9- 



-e- 




T 



V, 



Q • / 

Vi+l X i+1 



-e- 



T 



In the simplest case shown above, the (i + l)th y particle advances by the same 
amount as the ith x-particle. However in the event that this would result in it 
passing the position Xi, then it is blocked at that point, and the unused part of the 
increment (x^ — Xi) is passed to the (i + 2)th y particle as is shown beneath. 



Vi+l x i+ i 



Vi+2 



Xi+2 




T 



-e- 



L i+2 



We can also consider discrete analogues of these constructions. Set = 
A n n Z n /{2nN). Then, for r G {l/(2nN), 2/(2nN), ...,(N— l)/(2nN)} we may 
define a Markovian transition kernel on A^ as follows. Define 



1 if Z ~< r 



otherwise. 



and let 



(18) 



h N (x) 



n 



1 < / < m < n 
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where Xj = N N { _ n (xj + j /(2nN)). Then, as we shall see below, l^f, as a kernel on 
A% , admits h N as a strictly positive eigenfunction, and consequently 

(19) q r (x,X ) — foN J^j ( x i x )) 

defines a Markovian transition kernel, where 7^ > denotes the Perron-Frobenius 
eigenvalue of 1^ . We will call a Markov chain (X(k); k > 0) with these transition 
probabilities, an r- interlacing random walk on on . 

In the special case r = l/(2irN) the r-interlacing random walk on on A^ is 
closely related to non-coliding random walks on the circle, as considered in [25] for 
instance. Specifically if (X(k); k > 0) is the walk on A% then the process given by 

/ e iX- I (fc) )e iX 3 (fc) j _ _ ^ e iX n {k)\ 

where Xj(k) — N ^_ n (Xj (k) + j / (2nN)), gives a process of n non-coliding random 
walks on the circle { e ik/(2n(N+n) . k = o, 1, 2, . . . , iV + n - 1}. 

Fix r and s € {l/(2wN), 2/(2nN), . . . , (N - l)/(2nN)}. For u,v £ A% we may 
consider the two sets t^ v = r Ui „nZ n /(2nN) and t^' v = T^ v r\Z n /(2nN), where t u , v 
and t' u v are the polygons defined previously. The map <j> UiV carries l/ 1 /(2irN) into 
Z" /(2nN), and since we know it to be an isometry between t UjV and t' u v it must 
therefore restrict to a bijection between t^ v and t^' v . In particular the cardinalities 
of these two sets are equal. Since l^l^(u,v) = \t^ v \ = \t^' v \ = l^l^(u,v) we 
deduce that l$f and 1^ , and so and q^ 1 commute. This is useful in verifying 
our previous assertion that h is a eigenfunction of The case r = 1 may be 
easily deduced by direct calculation, see [25j for a similar calculation. Moreover 
the statespace A^ is irreducible for for r = 1 (though not in general), and thus 
the Perron-Frobenius eigenfunction of is unique up to multiplication by scalars. 
Consequently it follows from the commutation relation that h N is an eigenfunction 
of for every r. 

Using the fact that <j) uv is a bijection between t^ v and t^' v , we see that, 
substuting r„„ for r UjtI , for g r and so on, the proof of Proposition 13.11 applies 
verbatum in the case that X is an r interlacing random walk on A^ rather than 
A n . Similarly ip u<v is a bijection between t£ v and r^, and so Proposition [375] holds 
also in the discrete case. 
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We can extend the above constructions, by Kolmogorov consistency, to define 
a Markov chain {{X^ m '(k), m G Z}, k G Z) on the state space 

{x £ Al : a; (m) r<r z (m+1) , m G Z} 

with the following properties: 

(1) For each m G Z, (X( m )(fc),A: G Z) is an s-interlacing random walk 

(2) For each fc G Z, (X^{k),m G Z) is an r-interlacing random walk 

(3) For each m, & G Z, 

X< m >(* + 1) - 0x(™ +1 )( fc) ,x(™)( fc+ i)(^ (m) (fc))- 

In the above, we can replace by tp, and also consider the discrete versions. It is 
interesting to remark on the non-colliding random walk case, that is, the discrete 
model on A^ with r = s = 1/2ttN. In this case, the evolution of the above Markov 
chain is completely deterministic. Indeed, by ergodicity of the non-colliding random 
walk, for each k G Z, with probability one, there exist infinitely many m such that 
iW(ifc) g {((1 + c)/2irN, (n + c)/2%N), c G Z} and for these m there is 
only one allowable transition, namely to x!j m \k + 1) = X^ m \k) for j < n — 1 
and Xn 1l \k + 1) = X„(k) + 1/2-kN. By property (3) above, this determines 
X( m )(fc + 1) for all m G Z. 

4. Couplings of interlaced Brownian motions 

In this section we describe a Brownian motion construction that can considered 
as a scaling limit of the coupled random walks of the previous section. 

The Laplacian on A n with Dirichlet boundary conditions admits a unique non- 
negative eigenfunction, the function h defined previously at (fT2|) which corresponds 
to the greatest eigenvalue Ao = — n(n— l)(n+ 1)/12. A fe-Brownian motion on A n , 
is a diffusion with transition densities Qt given by 

(20) Q t (x,x')=e- x ° t ^Q t (x,x% 

where Q® are the transition densities of standard n-dimensional Brownian mo- 
tion killed on exiting A n , given explicitly by a continuous analogue of the Gessel- 
Zeilberger formula, |17|, 121] . If (X(k); k > 0) is an /i-Brownian motion in A n then 
^ e iXi(t)^ e iX 2 (t) ^ ^ e %x n (t)^ gj veg a p r0C ess of n non-coliding Brownian motions on 
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the circle, the same as arises as the eigenvalue process of Brownian motion in the 
unitary group, see |12j . 

We wish to construct a bivariate process (X, Y) with each of X and Y dis- 
tributed as /i-Brownian motions in A n , and with Y(t) ^ s X(t), for some fixed 
s G (0, 27r). The dynamics we have in mind is based on the map (f> of the previ- 
ous section. Y should be deterministically constructed from its starting point and 
the trajectory of X, with Yi(t) tracking Xi(t) except when this would cause the 
interlacing constraint to be broken. In the following we set 0o(t) = n (t). 

Proposition 4.1. Given a continuous path (x(t);t > 0) in A n , and a point 
y(0) € A n such that y(0) ^ s a;(0), there exist unique continuous paths (y(t);t > 0) 
in A n and (9(t); t > 0) in R" such that 

(i) Vi(t) =< s Xi(t) for allt> 0; 

(ii) Vi (t) - y<(0) = Xi(t) - Xi(0) + 6 l -i(t) - 9i(t); 

(iii) for i — 1,2, ...n, the real-valued process 9i(t) starts from 0i(0) = 0, is 
increasing, and the measure d9i(t) is supported on the set {t : yi + i(t) = 

Xi(t)}. 

Our principal tool in proving this proposition is another variant of the Skorohod 
reflection lemma. For I > 0, let A n {l) — {x £ R" : x± < x-i < . . . < x n < x\ + l\. 
We adapt our usual convention, letting = x\ + I and xo = x n — I. 

Lemma 4.2. Given a continuous path (u(i);t > 0) in R n starting from u(0) = 
0, and a point v(0) G A n (l) then there exists a unique pair of TV 1 -valued, continuous 
paths (v,9) starting from (v(0),0) with 

(i) v(t) £ A n (l) for every t>0, 

(ii) each component 9i (t) increasing with the measure d9i (t) supported on the 
set {t : Vi(t) = v i+ i{t)}, 

(iii) Vi(t) =w i (O)+u i (t)-0 i (i). 

Proof. Fix T > 0. The path u is uniformly continuous on [0, T], and so 
there exists e > such that J^. \ui(t) — Ui(s)\ < l/n 2 for all \s — 1\ < e. We prove 
existence and uniqueness over the time interval [0, e], and then repeat the argument 
over consecutive time intervals [e, 2e] etc. 



INTERLACED PROCESSES ON THE CIRCLE 19 

Observe that there must exist k E {1, 2, . . . , n} such that Vk+i(0) — Vk(0) > l/n. 
Set Vk(t) — Wfc(O) + Uk(t) and 9k(t) = for t € [0, e]. Then we the usual Skorohod 
lemma applied sucessively to give Vk-i,Vk-2> ■ ■ ■ j Vi, u n , . . . ffc+i, in particular we 
let 

8i(t) = sup(« 4+ i(*) - Ui(t) ~ Vi(0)y. 

s<t 

The choice of e is such that the Vk-i so constructed satisfies Vk(t) < Vk+\{t) for 
all t € [0,e], and consequently all the desired properties of v and 6* hold. For 
uniqueness, we note that any (v, 0) for which the desired properties holds must be 
equal to the one just constructed, which follows from the uniqueness for the usual 
Skorohod construction. □ 

Proof of Proposition 14.11 This is based on applying the Skorohod map- 
ping to the path u specified from x by Ui(t) = Xi(t) — Xi(Q). 

Choose v(0) G A n (2ir — s) so that Uj(0) — Uj_i(0) = t/i(0) — arj_i(0), and let 
(v, 6) be given by the Skorohod mapping for the domain A n (2Ti — s) with data u 
and v(0). Then for t > define 

Vi(t) = y 4 (0) + m(t) + - 0i(t). 

We claim that 

(i) t/(t) € A„ with :<„ x{t) for every i > 0, 

(ii) the measure d9i(t) is supported on the set {t : yi+\(t) — Xi(t)}. 

Calculate as follows. 

l«(t) - ar<_i(t) = (yi(0) + + - 0<(t)) - (x<_i(0) + Ui-i(t)) 

= ( Ui (0) + Ui(t) - 4 (*)) - («i_i(0) + Ui_i(t) - 0<_i(t)) = «i(t) - Wi_i(t). 

This is valid even if i = when we adhere to our conventions regarding vq etc. 
Assertion (ii) above follows since we see that {t : j/i(t) = Xi_i(t)} = {t : = 
Vj— and this latter set carries d9i—\. Also since Vi(t) > Vi—i(t) we obtain 
one part of the interlacing condition, namely, yi(t) > For the other part, 

consider the equality 



Xi (t) - yi{t) = a:<(0) - Vi(Q) + 0i(t) ~ 
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Since this quantity is initially Xi(0) — 2/i(0) > and decreases only when Q%-\ 
increases it follows that if there exists an instant t\ for which Xi{t\) — yi{t\) < 0, 
then there exists another instant, to, for which Xi(to) — yi(to) < and which 
belongs to the support of d0i-\. The latter implies that j/j(io) — ^i-i(io)) & n d thus 
Xi(to) < yi(to) — Xi_i(io) which contradicts x(to) £ A n . This proves existence. 
Uniqueness follows from the uniqueness statement in Proposition 14.21 □ 

By virtue of this result we may make the following definition. For y(0) G A n , 
let r a ( ) be the application which applied to an ^4„-valued path (x(t) > 0) returns 
the path (y(t);t > 0) specified by Propostion l4.ll The main result of this section 
is the following. 

Proposition 4.3. Let (X(t);t > 0) be an h-Brownian motion in A n , starting 
from a point X(0) having the distribution q s {y,dx) for some given y = y(0) € A n . 
Then the process Y = T y ^(X) is distributed as a h-Brownian motion in A n . 

The domain A n is unbounded, which is a nuisance when we come to the prob- 
ability, where life is easier if we have a finite invariant measure on the state space. 
Everything we do is invariant to shifts along the diagonal of R", and it is useful to 
project onto the hyperplane Ho = {x e R™ : x i = 0}- For x € R", the orthogo- 
nal projection of x onto Ho is given bymi-il where 1 € R" is the vector with 
every component equal to 1, and x — n^ 1 X{. If \X(f,)\ t > 0) is an /i-Brownian 
motion in A n , then its projection onto Ho is itself a diffusion process and indeed can 
be described as the /i-transform of an (n — l)-dimensional Brownian motion in Ho, 
killed on exiting A n n Ho- Introducing more generally H s = {x £ R" : Xi = s}, 
and the notion of an /i-Brownian motion on A n n H s , we have the following variant 
of Propostion 14.31 

Proposition 4.4. Let (X(t);t > 0) be an h-Brownian motion in A n n H s / 2 , 
starting from a point X(0) having the distribution q s {y,dx) for some given y — 
y(0) S A n H -H-s/2- Then the process Y = T y (n-\[X) is distributed as a h-Brownian 
motion in A n fl H_ s / 2 - 

Proposition ^. 3l is easily deduced from this variant using the fact that if (X(t);t > 
Oj is an /i-Brownian motion in A n , then the projection of X onto Hq is independent 
of the process n^ 1 ^ 2 Xi(t), which is a one-dimensional Brownian motion. 
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The results in the previous section were proved using the measure-preserving 
properties associated with the dynamics for a single update. Since here we are work- 
ing with continuous time processes such one time-step methods are not applicable. 
The use of the Skorohod lemma in the construction of T suggests that there may be 
a role to be played by a certain reflected Brownian motion. We describe next how, 
adapting the idea of Proposition 14. II slightly, we can construct interlaced processes 
X and Y from a reflected Brownian motion R in the domain HqC] A n (2 j n — s). Then 
it will turn out that time reversal properties of R can be used to prove Proposition 
PI 

Let E(s) — {(x, y) e (A n n H s / 2 ) x (A n n H_ s / 2 ) ■ y dis x}. We can introduce 
new co-ordinates on E(s) as follows. For (x, y) € E(s) we let f(x,y) be the unique 
(r,l) € H x H Q such that 

(21) r i+ i - n = y i+i - Xi, 

(22) -k = Xi+i - y i+ i, 

where r n+ \ = n + (27r — s) and l n +i = l\ + s. It is easily seen that / : E(s) — > 
(H a n A n (2ir - a)) x (H n A n (s)) is bijective. 

Now we construct a process in E(s) via these alternative co-ordinates. Begin by 
letting (U(t);t > 0) be a standard Brownian motion in R™ starting from zero, and 
let V(0) be an independent random variable, uniformly distributed in HqD A n (2-K — 
s). Let (V, &) be determined from U and V(0) by applying the Skorohod mapping 
for A n (2it — s) as given in Proposition ^. 21 Finally let R(t) be the projection of V(t) 
onto Ho, so R(t) — V(t) — V(t)l. The process (R(t);t > 0) is, by construction, a 
semimartingale reflecting Brownian motion in the polyhedron A n (2n — s) C\Hq. See 
[41j for the general theory of such processes. Next we introduce a process L also 
taking value in Ho, which is constructed out a random initial value L(0) together 
with the increasing processes 8^ for i — 1,2, ... ,n. Choose L(0) independent of R 
and uniformly distributed on A n (s) D Ho, and let L(t) be given by 

L(t) = Li(p)-nQ(t), 

where tt : R" — ► Ho denotes the projection onto Ho defined by irx = x — xl. Define 
the stopping time r = mi{t > : L(t) A n (s)}. Then for t < t, we define X(t) 
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andy(t) by (R(t),L(t)) = f(X(t),Y(t)) . The joint law of (X(tAr),Y(tAT);t > 0) 
may be described as follows. 

(1) (X(0),Y(0)) is uniformly distributed on E(s). 

(2) (X(tAr);t > 0) is distributed as a Brownian motion in Hq stopped at the 
instance it first leaves A ni and conditionally independent of F(0) given 
X(0). 

(3) y = ry ( o)X. 

We now turn to the time reversibility of R. For any vector x e R™ we will 
denote by the vector given by x\ — —x n -i + \. Note that if x e Hq n A n {2-K — s) 
then a;^ G J?o H A n (2n — s) also. We also define, for any x G R™, the vector via 
x\ = —x n -i — s/n for 1 = 1, 2, ... (n — 1) and x| = — a;„ + s(n — l)/n. Note that if 
x e H n A n (s) then i* e ff n A„(s) also. 

Proposition 4.5. Fix some constant T > 0. Let £/ie processes R and L be as 
above and let A be the event {L(t) e A n (s) for all t e [0, T]}. Then conditionally 
on A, 

(R(t),L(t);t e [0, T]) ^ (i?t(T - t), L*(T - t); t G [0, T]). 

PROOF. As remarked above i? is a scmimartingale reflected Brownian motion 
in the polyhedral domain A n (2n — s) fl Hq. Indeed it satisfies 



R(t) = R{0) + B{t) + -j=Y,e i {t)v i , 



where B is a standard Brownian motion in Hq and the vector v l describes the 
direction of reflection associated with the face Fi = {x G A n (2n — s) n i?o : = 
Let n J be the inward facing unit normal to this face. Then an easy calculation 
shows that v l , which is normalized so that the inner product n l ■ v l = 1, is given by 
v l = n' + q 1 where the jth component of the vector q l is given by 

V2/n-l/V2 \ij = i,i + l, 
\f2jn otherwise. 

We observe that the skew-symmetry condition, 



n l ■ q J + q 1 ■ n J = 0, 
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for all i ^ j, is met. Consequently by Theorem 1.2 of [41] . the reflected Brownian 
motion R is in duality relative to Lebesgue measure to another reflected Brownian 
motion on HonA n (2ir — s) with direction of reflection from the face Fi being n l — q l . 
It is not difficult to check that R) is such a reflected Brownian motion. Thus 

(R(t);t G [0,71) '= (R\T-t);te [0,T]). 

The process 2®i(t) is the local time of Ri+i(t)—Ri(t) at zero, and can be represented 
as 



1 /"* 

lim- / l(Ri +1 (s) - Ri(s) <e)ds. 



Now note that, 



/ l(i? l+ i( S )-i?,( S ) <e)ds'= / l(4 +1 (T- S )-4(T- S )<e)d S 
Jo Jo 

= / l[Rn^ i+ i(s) - Rn-i(s) < e)ds 

JT-t 



From this we deduce that the time reversal property extends to 

(R(t),e(t);te [0,T]) l = (R)(T-t),Q\T)-Q\T-t);t£ [Q,T]), 

where @\{t) = @ n -i for i = 1, 2, . . . , (n - 1) and 9^ = 6 n . 

Let F be a bounded path functional. Let v n (s) denote the Lebesgue measure 
of Hq H A n (s). Then using the time reversal property, 



E[F(R(t),L(t);te [0,T])1 A ] = 
1 



Vn(s)~ 



-E 



/ daF(R(t),a-irO(t);te [0, 7])l{a - 7r6(i) G A n (s);(e [0,T]} 
/ daF(R)(T-t) 8 J (T) +7r6 s (T-t);i G [0,T]) 

xlla-Tre^Tj + Tre^T-t) G 4(s);i£ [0,2*]} 



Wn(s) 
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Now we make the substitution a = a — 7r©"(T) to obtain 
1 



V n (s) 



E 



/ daF(R^(T-t),a + nQ i (T-t);te [0,T]) 

JHo 

xl{d + 7re«(T-<) e A n (s),te [0,T]} 



E[F(Jjt(T-t),i*(T-t);tG [0,T])1 A ] , 

where we have used = L*(0) + 7r6 s (t). □ 

Proof of Proposition 14.41 Let R, 9 and L be as above, and once again let 
A be the event that {L(t) G A n (s) for all t G [0,T]}. Recall the mapping / such 
that (R(t),L(t)) = f(X(t),Y(t)). It is easily verified that (yt(t), X+(t)) G B(s) 
and that 

(i2t(t),i*(t))=/(yt(t),jrt ( t)). 

Thus the preceeding time reversal result implies that, conditionally on A, 

(X(t),Y(t); t G [0, T}) ^ (yt( T _ t ), xt(T - G [0, T]). 

For the final step of the argument we consider X and 1" be as above and denote 
the governing measure by P. Then we let 

' oT -i A /i(y(o))/i(x(r))-p. 



% 

Under P, the equality in law, 

(X(t),Y(t);te[0,T}) ^ (Y^T-t),X\T-t);te[0,T]). 

holds unconditionally. Finally we note that under P, the distribution of (X(t),t £ 
[0,T]) and hence of (Y>(T—t), t G [0,T]) is that of a stationary /i-Brownian motion 
on A n n H s /2 But this latter law is invariant under time reversal, and its image 
under the conjugation x i— ► is the law of a stationary ft,-Brownian motion on 
A n n H_ s/2 . This is therefore the law of (Y(t),t G [0,T]). Conditioning on F(0) 
gives the statement of the proposition. □ 

In the case n = 2, the results of this section can be expressed in terms of 
Brownian motion in a compact interval. Let X — (X t ,t > 0) be a Brownian 
motion conditioned, in the sense of Doob, never to exit the interval [— p,p], where 
p > 0. Let y G [— p,p], a G [0, 2p] and suppose that the initial law of X is supported 
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in the interval [|y + a— p\ —p,p— \y+p — a\] with density proportional to cos(7ra/2p). 
Let Z be the image of the path X + (y — Xq + a)/2 under the Skorohod reflection 
map for the interval [0, a]. In other words, 



where L and U are the unique continuous, non-decreasing paths such that the points 
of increase of L occur only at times when Z t = 0, the points of increase of U occur 
only at times when Z t — a, and Z t S [0, a] for all t > 0. Then the process 



is a Brownian motion conditioned, in the sense of Doob, never to exit the interval 
[—p, p] . This is a special case of Proposition 14.41 Actually, in the statement of 
that proposition we have p = n, but this can be easily modified for general p. It 
is interesting to consider this statement when y = and p — > oo. Then X is a 
standard Brownian motion, initially uniformly distributed on the interval [—a, a]. 
The process Z is a reflected Brownian motion in [0, a] , initially uniformly distributed 
on [0, a]. The conclusion in this case is that Y is a standard Brownian motion 
started from zero. We remark that in this setting, if instead we take Xq = —a, 
then Y is a Brownian motion started from zero, conditioned (in an appropriate 
sense) to hit a before returning to zero. This is a straightforward consequence of 
the above result (for uniform initial law) and the fact (see [37] ) that, if we set 
T = mi{t > : Y t = a}, then the law of Xt is uniform on [—a, a]. Note that if we 
let a — > co in this case we recover Pitman's representation for the three-dimensional 
Bessel process. There are explicit formulae for the Skorohod reflection map for the 
compact interval [0,a] and hence for the process Y in the above discussion. Let 
fit) — Xt + (y — Xq + a)/2 and write f(s,t) — f(t) — f(s). A discrete version of the 
Skorohod problem was considered in [39 , from which we deduce the expressions 



Zt — max < sup min{/(r, t), a + inf f(s,t)}, min{/(t), a + inf f(s,t)} 



Z t =X t + (y-X + a)/2 + L t -U l 



Y t = y-X +Xt + 2(L t -Ut) 



t > 





An alternative formula was obtained in |26j . which yields 



Zt = 4>{t) - sup (4>(s) - a) A inf 4>{u) 

0<s<t [ s<u<t 
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where 

0(i) = f(t) + sup [-/(*)]+. 

0<s<t 

It could be interesting to relate the corresponding expressions for Yt to the Pitman 
transforms introduced in [4]. 

5. A bead model on the cylinder 

In this section it will be convenient to work with a slightly weaker notion of 
interlacing, defined as follows. For a, b € D n , write a -< b if 

ai < h < a% < ■■ • < a n < b n , 

and a >~ b if 

h < oi < &a < • • • < K < a n . 

For y = {e lai , . . . , e Mn } and z = {e lbl , . . . ,e lbn }, where a,b e £)„, write y -< 2 if 
either a -< & or a >- 6, and define 



Ej (^7 ~ a j ) + ^7T otherwise. 



Consider the Markov kernels defined, for q > 0, by 

(23) m,(|/,dz) = c" 1 /" |1 - e^r-Vpr^dsJdr, 

Jo 

where p r is defined by @ and 

f-271 



By Proposition ^. 1[ if we define 
then 



/ |1 -e ir r _1 9 r dr. 
Jo 



ify-<z, 

otherwise, 



(24) m q {y,dz)=i q 1 ^T^i (v,z)dz, 

where c g = c q /(n— 1)!. Recall that fj,(dx) = (2ir)~ n A(x) 2 dx is the probability 
measure on C n induced from Haar measure on U(n). The Markov chain with 
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transition density m q has fi as an invariant measure and, with respect to /z, has 
time-reversed transition probabilities 



m q (z,dy) = 5 « 1 ^fy X<? ^' : ;,// '" 



We can thus construct a two-sided stationary version of this Markov chain to obtain 
a probability measure a on C%, supported on configurations • • • -< x^ 1 -< x° -< 
x 1 < x 1 -< ■ ■ ■ . We will show that a defines a dctcrminantal point process on 
[0, 2tt) z . By stationarity it suffices to consider the restrictions a m to the cylinder 
sets C„ iTO := Cn 1 ' 2 '""'" 1 ^. Writing x — (x 1 , . . . ,x m ) and dx = dx 1 ■ ■ -dx m , 



(25) 



a m (dx) = fi(dx )m q (x , dx 2 ) ■ ■ ■ m q (x m 1 ,dx m ). 



Assume for the moment that q ^ 1. Define a function / : R — > C by 

^qgi(n—l)/2^u mod 2ir 
= 1 - (_l)n-l g 2 T • 

Lemma 5.1. for y = {e iai , . . . , e ia ™} and z = {e ibl , . . . ,e ibn }, where a,b £ D n , 

l q (y, z) = (l- (-I)"" V'Je' 2 ^ X^-W dct (/(6 fe - a,-))^- fc <„ • 
PROOF. Let c 7^ 1 and consider the n x n matrix W — {wjk) defined by 



w jk 



1 a,j < b k 
c a,j > bk- 



If a -< b, W consists of l's on and above the diagonal and c's below, so that 
det W = (1 — c)" -1 . If a y b, W consists of l's above the diagonal and c's on and 
below the diagonal, so that detW = c(c — l) n_1 . If neither a ~< b or a y b, then 
there must exist an index j such that, either dj = dj+i, or aj < dj+i < bk for all 
k, or bk < aj < dj + \ < bk+i for some k. In each of these cases, rows j and j + 1 of 
W are identical and hence det W = 0. Thus, 



(26) 



det W = { 



(l-c)"- 1 a^b 
c(c-l)™- 1 ayb 
otherwise. 
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Taking c = (ge^* -1 )/ 2 ) 2 * = (-1)"- V 7 ", we can write 



(1 - c)e 1 ^ det (f(bj - o fc )) 



l<j,k<n 



g Ej(6 3 -^)(l_ c )-(n-l) detW 

qUito-oi) a x b 
otherwise 



as required. 



□ 



For r = 1, ...,m— 1, define <?V r +i : [0, 27r) 2 — > C by r!r+ i(a,6) = f(b — a). 
Define <j> 0>1 : K x [0, 2tt) -> C and </> m , m+ i : [0,2tt) x I ^ C by 00,1 (a, 6) = e Mh 
and 4>m,m+i(a,b) = e~ lab . For r = l,...,m, write x r = {e Wl , . . . , e lo «}, where 
a r e D„, and set a° = a™ +1 = J — 1, for j = 1, . . . , n. 

Theorem 5.2. for 5 ^ 1, 



t (cte) = Zj Y[ det (</) r , r+ i(a^a; +1 )) 1< j k<n dx 



r=0 



where Z m = c™- x (l - (_i)«-i^)-(m-i) ( 27r )« 



PROOF. By l[2"I|) we can write 



,(«&) = a- (m - 1) (27r)-"A(a; 1 )A(a; m )X g (a; 1 ,a; 2 )---X (? (x m - 1 ,a; m )dx. 



Using the formula 



A(x 1 )A(a; w ) = det (e^-C^ 1 )/ 2 )^ det ( e -^-(™+D/ 2 K") 



l<j,/c<n 
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and Lemma EU we obtain 

a m (dx) = Z- 1 e i ^^M- a ^det(e i ^ n+1 ^ al A 

V / l<j,k<n 

m — 1 

x det ( e -<(WW +1 ) Yl Aet(^ r+1 (^ r a^)) l<j k<n dx 

= Z- 1 det f e i{j - 1)al A det f e^'- 1 ^ 1 ") 

\ ' l<3,k<n \ / l<j,k<n 

= Z- 1 [] ^t far+lW, < +1 )) 1 ^ k<n dx, 
r=0 

as required. □ 



m— 1 

X 
r=l 



COROLLARY 5.3. For any q > 0, the measure a defines a determinantal point 
process on [0, 2-7r) z with space-time correlation kernel given by 



K(r, a; s, b) 



hn^Ure^ r>s 

2tt 2^fcGZ\{0,...,n-l} 3fc e r < s 



where 



9k = \J " f(u)e- mk du) =i(k- ^ j - log?. 

Proof. For g ^ 1, this follows from Theorem 15.21 [22l Proposition 2.13] and 
a straightforward computation. The case q = 1 is obtained by continuity. □ 



Lemma ISTTl can be used to give a direct proof of (f24|) . and hence Proposition ^. II 
PROOF of Proposition 12.11 The characters x\ are given, for y = {e lai , . . . , e w ™}, 

by 

Xa(») = r(f A(y)" 1 det (e^ Qfc ) > 
where /i = X + p and 

( n— In — 1 n — 1 n — 1 



2 ' 2 2 '2 

Using Lemma 15. II and the Cauchy-Binct formula, we obtain 

/ ^^(»^)xA(«)^ = (i-(-i) w - 1 s to )IlH^-iog?)~ 1 XA(i/). 
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On the other hand, writing x r — {e lr , 1,1,..., 1}, an easy calculation shows that 

f V |1 _ eir\ q rXXp± dr = {n __ _ (.l)""!^) JJ (-i _ lQgq) -X 

Jo d x Y 

and so, by J7]), 

m qX x = c-\l - (-lr-y*) - log?)- x XA- 

j 

Since {xa,^ G f2„} is a basis for L 2 {C ni n), this implies □ 

Analogous results to those presented in this section can be obtained for the 
discrete version of this model, which is equivalent to considering a certain family 
of Gibbs measures on rhombic tilings of the cylinder. For more details, see [20] . 
The couplings defined in Section 3 are quite useful in this setting, where the group- 
theoretic considerations of Section 2 no longer apply. For example, they can be 
used to prove that the discrete analogues of the interlacing operators {T q , q > 0} 
commute with each other. In this setting, the symmetric functions 

(qi,...,q k ) H->J gi ■■■I qk (y,z) 

are essentially the cylindrical skew Schur functions discussed in the papers [27\ 136) . 

Finally, we remark that, in the case q = 1, the probability measure a 2 defined 
by (|25p also arises naturally in random matrix theory. The probability measures on 
C„ given by fi(dx) — (2ir)~ n A(x) 2 dx and A~ 1 A(x)dx, where A n is a normalisation 
constant, are known, respectively, as the circular unitary ensemble and circular 
orthogonal ensemble. It is a classical result, which was conjectured by Dyson [11] 
and subsequently proved by Gunson 20J , that the set of alternate eigenvalues from a 
superposition of two independent draws from the circular orthogonal ensemble, are 
distributed according to the circular unitary ensemble. Moreover, the joint law of 
the 'even' and 'odd' eigenvalues has probability density on C n x C n proportional to 
A(y)A(z)2"i(y, z), which is the same as a% the joint distribution at two consecutive 
times of the stationary Markov chain with transition kernel m q (defined by (|2"3"]) ) 
in the case q=l. 
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